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Abstract 

After the work of Rindler and Ishak, it is now well established that the 
bending of light is influenced by the cosmological constant A appearing in the 
Schwarzschild-de Sitter spacetime. We show that their method, when applied 
to the galactic halo gravity parametrized by a constant 7, yields exactly the 
same 7— correction to Schwarzschild bending as obtained by standard methods. 
Different cases are analyzed, which include some corrections to the special cases 
considered in the original paper by Rindler and Ishak. 



I. Introduction 

It has long been believed that light bending in the Schwarzschild-de Sitter 
(SdS) spacetime is uninfluenced by the cosmological constant A appearing in the 
metric. The reason is that A cancels out of the second order null geodesic equa- 
tion — naturally, the light trajectory too does not contain it. On the other hand, 
A appears in the first order differential equation, only its further differentiation 
removes the constant A from the second order differential equation. Obviously, 
for the sake of consistency, the perturbative solution of second order equation 
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must satisfy the first order equation as well, which would then yield a relation 
among the involved constants, one of which is the impact parameter b. This 
would in turm imply that the removed A will reappear in the geometric light 
trajectory, hence in the light bending as well. We shall address this question 
in more detail elsewhere. Here we shall limit ourselves only to the perturbative 
solution of the second order equation expressed customarily in terms of closest 
approach distance 

Recently, a new method has been proposed by Rindler and Ishak [1] that 
combines the standard perturbative solution with an invariant geometric defini- 
tion of the bending angle. The method reveals that the Schwarzschild bending 
caused by a lens (say, a galactic cluster of mass M) is decreased by repulsive 
bending due to A. This work has instantly generated a lot of interest among 
the gravity physics community (see, e.g., [3,4,5,6,7,8,9,10]). It can be fairly said 
that the original work by Rindler and Ishak (as well as its extension [2] to the 
Einstein- Strauss vacuole model) has established beyond doubt that there is a 
A— dependent effect on light bending contrary to previous belief, although its 
interpretation in the cosmological set up is still open to further discussion and 
research (see the recent review [11]). In view of this new wisdom, it would be 
interesting to apply the Rindler-Ishak method to the domain of galactic halo 
gravity. For this purpose, an excellent example seems to be the Mannhcim- 
Kazanas-de Sitter (MKdS) solution of Weyl gravity with a conformal parameter 
7, interpreted as a player in the halo gravity. The Weyl theory is based on the 
15-parameter group of conformal invariance and it attempts to resolve the dark 
matter/dark energy problem without hypothesizing them. It should be noted 
however that the MKdS solution has its own merits and demerits, but their 
discussion is beyond the scope of this paper. Our interest lies only in finding 
if the Rindler-Ishak method can reproduce the first order effect of 7 on light 
bending already known in the literature. We show that it indeed does, which 
can be treated as yet another success of the method. 

The purpose of the present paper is to implement in a more general spacetime 
the Rindler-Ishak method to one higher order in M than considered originally, 
and calculate the effect of 7 on the bending of light rays. It turns out that 
the method delivers the exact first order 7— term in addition to revealing new 
interplays among the constants A, M and 7. Different cases are discussed and 
some needed corrections in [1] are pointed out. 

The contents are organized as follows: In Sec. II, we derive the geodesic 
equation in the MKdS solution. In Sec. Ill, we work out the bending of light 
rays following the Rindler-Ishak method. We discuss specific cases in Sec. IV 
while Sec.V summarizes the paper. There are three appendices. 

1 Following Weinberg [16], we shall take in the light orbit the closest approach distance R in 
preference to the impact parameter b. This then implies that we can carry on with the second 
order differential equation with its perturbative solution involving R. It should be mentioned 
that Weinberg's integral too eventually involves only R (not b) in which A disappears (see 
[13]). The method has been applied only to the asymptotically flat spacetimes, which is not 
the case here. Thus, a new method, such as that of Rindler and Ishak, seems more appropriate, 
which we follow here. 
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II. Geodesic equation 

An interesting solution of the Weyl gravity field equations is the MKdS 
metric given by [12,13] (in units 8nG = c = 1): 

dr 2 = B(r)dt 2 1 —dr 2 -r 2 (de 2 +sm 2 0d<p 2 ), B(r) = 1- — + 7 r-fcr 2 , (1) 

B(r) r 

where M is the luminous central mass, k and 7 are constants. The numerical 
value of k — A/3 = 0.43 x 10~ 56 cm~ 2 , and 7 is of the order of inverse Hubble 
length. However, we caution that there is a reported ambiguity both in the 
magnitude and sign of 7. By analyzing the fiat rotation curve data, Mannheim 
and Kazanas [12] fix it to be positive, 7 s=s + 1CP 28 cm -1 , while Pireaux [14] 
argues for 7 ps — 10~ 33 cm -1 . Edcry and Paranjape [13] conclude from the time 
delay data in the halo that 7 ~ — 7 x 10~ 28 cm -1 . We emphasize that we are 
free to adduce any sign to 7 in the sequel but for definitcness, we choose 7 > in 
the present work. Such choice is neither mandatory nor essential for the present 
work. 

Denoting u = 1/r, we derive the following path equation for a test particle 
of mass mo on the equatorial plane 6 — ir/2 as follows: 

d2u „„„ 2 7 M 1 / 2fc\ 

where h = the angular momentum per unit test mass. For photon, m = 
=> h — > 00 and one ends up with the null geodesic equation without k making 
its appearance: 

d 2 u „ , , 9 7 

— = -u + 3Mu 2 - 1 
dtp z 2 

Here we find that a cancellation of k similar to that in the SdS case (7 = 0) 
occurs despite the presence of a nonzero 7 in the metric. Exactly like in the 
SdS case, one would now tend to believe that the bending of light in the MKdS 
case would not be the influenced by k to any order but this is not the case. 

III. Rindler-Ishak method 

We shall follow the three principles adopted by Rindler and Ishak [1]: (i) 
The method is originally applied to the SdS metric to show that, despite the 
non-appearance of k in the null geodesic, its effect still appears in light bending. 
Hence to preserve the essence of the method, we shall retain k 7^ except in 
special cases, (ii) With M 7^ 0, k 7^ in the metric, the limit r — > 00 makes no 
sense. The only intrinsically characterized r value replacing r — > 00 is the one 
at ip — 0. The measurable quantities are the various ip angles that the photon 
orbit makes with successive coordinate planes ip — const. We shall qualitatively 
verify the results for ip 7^ as well, say, at ip — ir/4. (iii) We shall use the 
perturbative solution for - up to order M 2 and the resulting deflection angle ip. 

Let us develop the basic equations now. Although the MKdS metric is more 
general than the SdS metric, we shall show that the influence of 7 appears in 



^ = -u + 3Mu 2 -l (3) 
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the light bending together with that of k including terms mixing up the two. 
The light path equation (3) in the zeroth order is 



d U ° , T (A\ 

1^ +Ua= 2 (4) 

and its exact solution is 

COS (£> 7 

«o = -^-? (5) 

where i? is the distance of closest approach to the origin (where the lens is). 
This solution is to be used as the zeroth order approximation. Following the 
usual method of small perturbations [15], we want to derive the solution as 

u = u Q + m + u 2 (6) 

where u\ and u 2 respectively satisfy 

d 2 ui 



dip 2 

d 2 u 2 
~d^ 



+ m = 3Mw 2 (7) 



+ u 2 = 6Mu ui. (8) 



The exact solutions are 

M 
4R2 



u\ = [6 + 3R 2 "/ 2 — 6-R7 cos tp — 2 cos 2ip — 6Rjtp sin ip\ (9) 



3M 



2 



U 2 



• [{ 10 + 3i? 2 7 2 (5 - 2ip 2 ) } cos Lp - 12,R 7 (4 + i? 2 7 2 ) 



16i? 3 ' 

+ cos3(/3 + 20<y9sin<y9 + 16i?7Cos2<^ + 30R 2 ^y 2 tpsmtp + SRjtp sin 2<p] . (10) 



Formally changing ip — > | — i^, the final solution up to second order in M can 
be written as 

1 sin ip 7 
?i =r = ^T-2 

+ ^ [ 6 + 3i? V _ 3i?7(7r - 2ip) cos ip + 2 cos 2^ - 6.R7 sin 93] 

3M 2 
~ 32i? 3 

- 20 sin ip - 30i? 2 7 2 sin tp + 3tt 2 R 2 j 2 sin tp - 12ti\R W sin (p 

+ 12i? 2 7 V sin tp - 87ri?7 sin 2ip + WRjtp sin 2ip + 2 sin 3<p] . (11) 

This pcrturbative expansion holds only for small u or large r. Thus we are 
considering parameters M, R and 7 such that ^ < 1 and 7 i? < 1. When 
7 = 0, it may be verified that one recovers the equation for light trajectory up 
to M 2 order in the Schwarzschild metric. The following clarification be noted: 



[96i?7 + 24i?V - 10(2 + 3R 2 Y)(ir - 2ip) cos tp + 32Rjcos2tp 
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In the expression for u in Ref. [2], a trivial term 10 simp is deleted while another 
trivial term 107rcos<^ is retained. From the standpoint of generality, we keep 
our exact orbit equation as it is, but show that the presence or absence of such 
trivial terms do not lead to any difference whatsoever in the final result (See 
Appendix C). 

The method of Rindler and Ishak is based on the invariant formula for the 
cosine of the angle ip between two coordinate directions d and S such that 

C0S ^ = j^W^m^' (12) 

Differentiating Eq.(ll) with respect to ip, and denoting ^ = A(r,<p), we get 

A(r,ip) = (-r 2 )x 

[g^g(32i? 2 - 3M 2 {20 + 3i?V(10 + (tt - 2<^) 2 } 

+ 32Mi?(9M 7 - 2) sin ip) - 6M{3M cos 3p - 8Mi? 7 (7r - 2ip) cos 2ip 
+ (10M - 4i? 2 7 + 9Mi? 2 7 2 )(7r - 2ip) sinp}}. (13) 

Eq.(12) then yields [1] 

CO sV= — r - (14) 

[A 2 + B(r)r 2 ]i 

In a more convenient form, the final Rindler-Ishak expression for ip to be used 
is [see (iii) above] 

B i/2 r 

tan ^ =- j]4p (15) 

The bending angle is defined by e = ip — ip. The basic ingredients are Eqs.(l), 
(11), (13) to be plugged into Eq.(15). We shall now discuss specific cases. 
IV. Specific cases 

Case 1: ip — 0, M ^ 0, k ^ [This is Rindler-Ishak choice; see also (ii)]. 
We get from the orbit Eq.(ll) 

r = 16R 3 /[4MR(8 - 37ri? 7 + 3i? 2 7 2 ) - 8i? 3 7 
+ 3M 2 {5tt(2 + 3i? 2 7 2 ) - 4i? 7 (16 + 3i? 2 7 2 )}]. (16) 

The SdS value of r in the first order (in M) solution in Ref.[l] is exactly 
which can be recovered by neglecting the M 2 term in Eq.(16) and setting 7 = 0. 
For the unbound orbits associated with lensing the distance of closest approach 
of a light wave to a galaxy will be further from the center of the galaxy than 
the matter orbiting in it. Consequently, for practical purposes for lensing one 
should consider the halo regime where jR > ^ holds. We have been considering 
expansion (11) in the Schwarzschild weak field regime due to M so that ^ < 1. 
We also assume 7 i? < 1 so that (jR) 2 terms can be neglected. 

We shall implement these approximations on the results derived from the or- 
bit Eq.(ll). Let us identify the leading order term in r taking into account prac- 
tical data, say, for a typical galactic cluster Abell 2744, for which the accurately 
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observed central mass and Einstein radius respectively are M = 2.90 x 10 18 cm 
and i?E = 2.97 x 10 23 cm [2]. The function r in Eq.(16) expands term by term 
to first order in 7 as 

R 2 i? 4 7 157ri? 37Ti? 3 7 2 457r 2 i? 2 7 

r = 1 + 3R -7 - 

2M 8M 2 32 64M ' 2048 

225tt 2 M 45irMR~/ 6757r 3 Mi? 7 ^.,, 2n 

+ ^ r J1 + —m J1+0 ^ (17) 



For light moving in the halo regime, R > Re, say, R = 3 x 10 23 cm and with 
7 = 7 x 10~ 28 cm -1 [12], it can be verified that the condition jR > ^ holds 

and that the leading value of r is ■ Putting this value of into B(r), with 

k = 0.43 x 10~ 56 cm~ 2 [12], we get the signature protection 

B (S) >0 ' <18) 

Putting the value of r from Eq.(16) into Eq.(13), we get 

A = 8i? 3 [M 2 (78 - 48^7 + 90i? 2 7 2 + 9ir 2 R 2 j 2 ) - 32i? 2 ]/[4Mi?(8 - SttRj 
+ 3i? 2 7 2 ) - 8i? 3 7 + 3M 2 {5tt(2 + 3i? 2 7 2 ) - 4^(16 + 3i? 2 7 2 )}] 2 . (19) 

Ignoring M 2 terms and putting 7 = above, we recover the value of \A\ in the 
SdS spacetimc. The one sided bending angle is given by e = ip — tp and let us 
calculate e — tp — ^l> when ip = 0. Putting in Eq.(15) the values for B(r) from 
Eq.(l), r from Eq.(16) and A from Eq.(19), we get 



tanr/i = y\j 1 — + ^ r — kr 2 j x 

2[AMR{S - 3wR-y + 3i? 2 7 2 ) - 8i? 3 7 + 3M 2 {5tt(2 + 3i? 2 7 2 ) 

- 4i?7(16 + 3i? 2 7 2 )}]/[M 2 (78 - 487ri? 7 + 90i? 2 7 2 + 9ir 2 R 2 j 2 ) - 32R 2 }. 

(20) 

This is the exact formula for light deflection but is rather unilluminating. There- 
fore, successively expanding in the first order of 7, and in the second order in 
M, we find, for a small angle ipo (or, tani/'o — ipo), the following expression (see 
Appendix A for details, see also Appendix C): 



2M 



, WirM kR A 
1 + ■ 



16R 8M 2 



R 3ttM 455M 2 4M 
+ —r- + ^ + 



2 4 32i? kR 3 



(21) 



where we have retained only the leading order terms in the coefficient of 2M. 
and 7 assuming ^ < 1. We find that all terms in the last third bracket are 
positive, meaning that the effect of 7 is to diminish the Schwarzschild bending 
even up to second order in M. This is a core result of this paper. 
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We find that Eq.(21) reproduces the exact term — ^ obtained earlier in the 
literature by Edery and Paranjape [13] using Weinberg's method. Thus, it is 
clear that the Rindler-Ishak method is not only correct but also tells us more 
in the form of other new terms showing, in the last bracket, a clear interplay of 
7, M and k, each bearing its own physical meaning. When the light orbits in 
the halo region, R> Re, it was already computed that the 7— term dominates 
over the Schwarzschild term: ^ > 2|r-, i.e., the parameter 7 plays the main 
role in the halo. At R > Re, the effect of M ^ would be practically negligible 
so that only repulsive bending — ^ will occur. However, the (7/fc)— dependent 
term could be comparable to the pure 7— term in spite of the fact that in the 
halo region, << 1. We need realistic values to see that fact. For the values 

stated around Eq.(17), we get: ^f- ~ 1(T 4 , ^£f- ~ 1CT 4 which shows that 
the pure 7— effect is of comparable magnitude with the combined effect of 7, 
M and k. The other terms are: ^ ~ 1CT 5 and ~ lCT 7 and are at least 

an order of magnitude less than ^-j 2 -. So, the Schwarzschild bending seems 
to be dwarfed by the 7— terms. However, as we see, this conclusion is sensitive 
to the exact magnitude of 7. The problem is that, while the values of M, Re, 
and k are observationally known with adequate accuracy [2] , the exact sign and 
value of 7 are not yet conclusively known. 

Case 2: <p ^ 0, M = 0, k 7^ 0. The exact solution of the null geodesic is 
- = ^j^- — 3r, which for small angle ip reduces to 

l - = ^-l. (22) 
r R 2 V ' 

Then the angle ip corresponding to r = r max is 

7.R kR 

<^(r max ) = V + — . (23) 
2 7 

According to the present method [see (hi)], we have, differentiating (22), 

^(w) = (-rLx)^=-^2- (24) 

Putting the values of r max = ^, \A(r max )\ and B = B(r max ) = 1 in (15) for 
small ip, we obtain, 

kR 

^max) = • (25) 

7 

The deflection is given by 

7.R 

£ = ^(^max) - ^(^max) = ~ • (26) 

The final result (26) is independent of k. Because of this independence, the 
deflection in the pure conformal gravity with metric potential B(r) = 1 + "/r 
also gives the same deflection, viz., e — — ^ (See Appendix B for an alternative 
derivation). 
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Case 3: p = 7r/4 [Rindler-Ishak choice], M 7^ 0, k 7^ 0. Below we make 
some corrections. The Rindler-Ishak formula gives the bending angle e as 

. , . tan ip — tan p tan ip — 1 , . 

e = ip - ip n tim(ip - ip) = ^— — = - - -. (27) 

1 -I- tan ip tan p l+tan-0 

Proceeding in the similar manner as above, we can find the value of tan ip and 
thence of e. Surprisingly the coefficient of becomes identically zero. We 
do not give detailed expressions here but only state the final result to leading 
order: 

V2M 3irjM 7 2 i? 2 fc 7 2 i? 4 kR 2 

£=_ R 8 8 4 2~' (28) 

Once again the deflection due to 7 is negative. Incidentally, the first term is 
slightly different from that of Rindler and Ishak [their Eq.(19)]. Also, for 7 = 0, 
their expression [Eq.(18)J for tan?/> slightly modifies to 

2V2M , , 

tan^ = l + -^ kR 2 . (29) 

R 

These modifications by no means alter their demonstration of the repulsive 
effect of k (in fact it is exactly the same). Even the pattern (loosely speaking) 
of "division by 2" in the last two pieces of taxnp results in the corresponding 
pieces in e, just as it is in Ref. [1]. The minor changes in the coefficient of M 
have come about because we have used the path equation at p = 7r/4 to obtain 
r = V2R - 3M + 0{M 2 ) rather than r = V2R. 
Case 4: p = tt/4, M = 0, k ^ 0. Then 

2R t 2V2R , , / 7 2 i? 2 + 4fci? 2 - 2 \ 1/2 

r = —= , A = — = , tan p = 1, tanw = — = 

V2- 7 i? (V2- 7 i?) 2 \4V2jR-2j 2 R 2 -4/ 

(30) 

For small deflections, expanding successively in powers of 7 and first power of 
k, we get 

tan ip- tan p j 2 R 2 fc 7 2 i? 4 /ci? 2 ,„ , 

e = ~ . 31 

1 + tan ip tan p 8 4 2 v ' 

This result can also be obtained directly from Eq.(28) by putting M = 0. 
V. Summary 

As shown in Eq.(3), in the MKdS gravity too, which is more general than the 
SdS gravity, the constant k cancels out of the light orbit equation even though 
7 7^ 0, the latter fact distinguishing the MKdS metric from the SdS metric. We 
should also remember that the corresponding parent theories are generically 
very different; one is fourth order and the other is the second order gravity. 
Nevertheless, in view of similar cancellation in the two metrics, we investigated 
the applicability of the Rindler-Ishak method taking the MKdS solution as an 
example. This is a nice example because Weyl conformal gravity accommodates 
the successes of Schwarzschild gravity and has been the subject of active research 
for several years. We should re-emphasize that we are exclusively concerned here 
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with the efficacy of the Rindler-Ishak method and not with the well discussed 
but inconclusive values of and other difficulties associated with the 7— term. 

We first derived the exact null geodesic equation in the MKdS gravity. Next, 
we perturbatively solved the equation up to the order M 2 though the zeroth or- 
der equation is different from that in standard Schwarzschild gravity. Whatever 
follows below are the results obtained after the detailed solution is plugged into 
the Rindler-Ishak procedure, which we have faithfully implemented. 

We should note that, generally speaking, none of the quantities k, 7 and M 
should be zero. (We can nonetheless set one or the other to zero as limiting 
cases). Then Eq.(21) at once shows the influence of k both in the Schwarzschild 
and conformal sector. The equation nicely reproduces the attractive Schwarzschild 
bending due to M > as well as the repulsive bending due to cosmological con- 
stant k > and the Weyl conformal parameter 7 > 0. Eq.(21) combines in 
one place deflections at various orders obtained by Bodenner & Will [15], by 
Rindler & Ishak [1] and by Edery & Paranjape [13]. Because of M 2 order, there 
appeared new terms in the conformal sector. In particular, there is a second 
order correction — 3lT J M = — 2zi x jR x ^ to the first order term — ^ obtained 
in Ref.[13], which also diminishes the Schwarzschild one way bending. 

There is absolutely no problem in the case ip — tt/4, k — 0, M — [Eq.(28)] 
because the expression for e does not blow off. In this case, the term — ^- of 

course does not appear, but one notices that the effect is — 7 8 , which remains 
unaltered by the choice of sign for 7. In either case (tp = or 7r/4), the conformal 
parameter 7 can always be set to zero at will, but then one ends up with the 
already discussed SdS case. For the case M = 0, we determine the maximum 
allowed value of r = r max and correspondingly determine e = — ^ [Eq.(26)]. 
A direct integration given in Appendix B also supports this result. Since e is 
independent of k, this then is also the deflection in the pure conformal gravity 
defined by M — 0, k — so that B(r) = 1 + jr. Overall, the conclusion is 
that the Rindler-Ishak method can be applied to more general situation than 
SdS and that it leads to the same result up to second order as obtained by 
conventional perturbative method (See Appendix C). 

The implication of the last term in Eq.(21), viz., — -7^753-, is very interesting 
and seems to provide the background for a certain postulate. To have some idea 
how, let us estimate the terms in Eq.(21) for a light ray grazing a stellar sized 
massive object, say, the Sun itself, and that 7 w 10~ 28 cm~ 1 , k w 10~ 56 cm~ 2 , 
the values being already independently estimated. We then have the following 
numerology: 

Mq = 1.48 x 10 5 cm, R Q = 6.96 x 10 10 cm, (32) 



so that 

2Mf 



i£ = 424 x 1(r6 ' ^kf - 2 - 65 x 10 ' - im! - - 10 

2^ « -10-" - „ _!„-» _ ^ „ _ 10 . 

2 '4 ' kR% 



(33) 
(34) 
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We find that, while all the other terms above are quite small compared to -j^-, 
the spoiling term is the very last one. If we had limited ourselves to the direct 
integration of the null trajectory in which k does not appear, we would have 
missed this term. This term dominates at the limb of the Sun, giving rise to 
a total bending, ipo m — 10 6 , which is negative and means repulsion 10 12 times 
bigger than the first order Schwarzschild attraction! Certainly this is contrary 
to our experience. Moreover, with M = Mq, for any R > Rq, one finds that 
the values of ipo continue to remain only negative. Without that spoiling term, 
however, the effect of 7 would indeed be negligible near the Sun so that ipo 
would not appreciably differ from the positive Schwarzschild value, which docs 
not appear to be the case here. 

One possible viewpoint is to postulate that, at the solar scale, MKdS is 
predominantly only SdS [negligible 7— effect like in the first two terms in (34)], 
while its genuine applicability lies at the galactic scale. In fact, the deflection 
on the galactic cluster scale tells a quite different story: Consider again Abell 
2744 [2], with £ w 10 28 cm, the term -^jf- « -1(T 4 , which is comparable 

in order of magnitude with that of — ^ and is at most an order of magnitude 
higher than the attractive Schwarzschild term if we believe in the value of 7 used 
here. In any case, the overall repulsion thus obtained in the halo can always be 
converted to the desired attractive bending if the numerical value of 7 is slightly 
altered for which possibilities certainly exist. The main conclusion is that, with 

7 ~ in the solar system, the huge repulsion (— 4 ^a ~ — 10 6 ) term can be 

avoided. This implies that 7 could be physically relevant only on the galactic 
cluster scale and not on the solar scale. The fact that 7 operates in the galactic 
halo has been conjectured in the literature, but here we find its support from a 
completely different viewpoint, viz., from the Rindler-Ishak bending. 
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Appendix A 

For the benefit of the readers we give below the exact steps leading to the 
expression (21). After putting the value of r from Eq.(16) into Eq.(20) and 
simplifying, we obtain 



tan^o 



1 

S 



2(P + Q){l + 7 M 



251 i A- 



8AF 



167iT 

P+Q (P+Q) 



2W 

1/2 



A? 6 \ 



where 



P = 4MR{8 - 37r 7 i? + 3 7 2 i? 2 ) - 8 7 i? 3 
Q = 3M 2 {5tt(2 + 3 7 2 i? 2 ) - 4 7 i?(16 + 3 7 2 i? 2 )} 
S = M 2 (78 - 487ri? 7 + 90i? 2 7 2 + 97r 2 i? 2 7 2 ) - 32i? 2 . 



(Al) 

(A2) 
(A3) 
(A4) 



In order to extract the contribution due purely to 7 , we expand the right hand 
side in the first power of 7 , and get 



tan -00 



2TU 



78M 2 - 32R 2 



2 7 T 



M 



3ttM 2 24A/ 3 16R 3 



265kR 7 W 



2R 



R 2 



T M 3 (15ttM + 16R) 3 



(A5) 



where 



U= 1- 



T = ZOirAr + 32MR 
4M 2 15ttM 3 256/ci? 6xl/2 



v = u 



156 M 2 - 64i? 2 - 



R 2 4i? 3 T 2 

2RW 247ri?M 3 (157rM + 16R) 



39 M 2 - 16i? 2 (39Af 2 - 16i? 2 ) 2 

W ee 48M 2 + SttMR + 2R 2 . 



(A6) 
(A7) 
(A8) 
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We first expand the second term in (A5) in powers of M obtaining 

455M 2 4M 2 " 1 



MKdS 



"7 



R 3irM 
2 + — 



(A9) 



32R kR 3 

Next let us put 7 = in (A5) in order to obtain pure Schwarzschild terms from 

2TU 



tan^ dS 



78M 2 - 32i? 2 ' 



(A10) 



The deflection is already known to be positive. To ensure it, Rindlcr and Ishak 
in their SdS treatment prescribed \A\ instead of just A. This prescription is the 
same as changing the denominator of (A10) into 32i? 2 — 78M 2 because R » M 
and TU > 0. For small Vo, and with U ~ 1 - - 1^1 - the first 

term in (A5) then results in 

2(30ttM 2 + 32Mi?) 
^SdS = 



1 _ 2M' 
1 TF" 



128fcft° 



32i? 2 - 78M 2 
which, when expanded in powers of M, yields 

15ttM kR 4 



,SdS 2M 
^ ~ IT 



1 + 



16i? 8M 2 



The total one way deflection is of course 

^ = ^SdS + V ,MKdS 

which is just the Eq.(21) in the text. 

Appendix B 

Consider the null geodesic equation when M = 0: 



dip' 



7 

- U -2 



which has an exact solution 



The metric for M — is 



1 sin ip 7 
r = ~R~ ~ 2' 



B(r) = 1 + 7r — fcr 2 . 
Weinberg's method allows a direct integration giving the deflection as 



f 



r i {l+ 1 R) 2 3 

- r z - 7r d 



R? 



-1/2 



r 

i? 2 



'1/2 



7i? 



-1/2 



rdr 
R + r' 



(All) 



(A12) 



(A13) 



(Bl) 

(B2) 
(B3) 



(B4) 
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The last line is obtained after expanding the integrand in small jR. As noted 
by Edery and Paranjape [13], k has cancelled out of the integrand. It is to be 
expected since the path equation does not contain k either. Therefore, even for 
pure conformal gravity, when B(r) = 1 + -yr, the above result holds true. The 
integration in (B4) yields 

A,-f=-f. (B5) 

This result supports the fact that pure conformal gravity is repulsive. The factor 
§ appears due to the conventional definition of the bending angle. This factor 
does not appear according to Rindlcr-Ishak definition [1] of the azimuthal angle, 
but the final result is always the same once it is corrected for, which is indeed 
the case as shown in Eq.(26) in the text. 

Appendix C 

It is well known that the structure of complete solution of an ordinary dif- 
ferential solution has two parts: The Characteristic Function (CF) is one part 
when the right hand side is zero and the Particular Integral (PI) is the other 
part when the same is nonzero. One might argue that terms proportional to 
simp, cos ip of the perturbativc differential equation (8) [by the same token, also 
of Eq.(7)] are trivial CFs and hence should be discarded from the nontrivial PI 

part of the solutions. Indeed, in the SdS light trajectory (with 7 = 0), the full 

.2 

solution of the perturbed zeroth order equation + u = contains only 
u cc sin ip, cos tp as CFs and no PI because the right hand side is identically 
zero. So one might wish to avoid their repetitive occurrence in the other two 
higher order perturbed solutions. This can be easily achieved by choosing the 
constants of integration. However, in the case 7 / 0, the full zeroth order 
solution is 

u = CF+PI = (Ci sin ip + C 2 cos tp) + - I) . (CI) 

We can choose C\ — so that the reduced CF is C2 cos tp. If we do not want 
its repetition in the nontrivial PI part (^jp — 77), there are two ways: We can 
simply delete co ^ p from it and choose C2 = 1 / R assuming R to be the distance 
of closest approach. Alternatively, we can set all the arbitrary CF constants to 
zero and retain only the PI as it is, even if it includes the trivial term. In either 
way, we would end up with the same uo — — ^. We want to examine the 
conventional Bodenner-Will method [15] for the MKdS solution using both the 
ways (viz., with or without CFs in PI). 

The structure of the exact full solution (CF+PI) of the three second order 
perturbative equations (4), (7) and (8) would involve six arbitrary integration 
constants C\, C2, C3, C4, C5 and Ce so that in the Schwarzschild case, we find 

u = [Ci sin f + C 2 cos tp 
+ (C 3 + Ci) sin ip + (C 5 + C 6 ) cos tp] 
15M 2 M 3M 2 
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The form of the solution considered by Bodenner and Will [9] is equivalent to 
choosing 

1 1 'iM 2 

C 1 =0,C 2 = -,C 3 + C 4 ^0,C i + C 6 = -^ r , (C3) 

which removes the trivial CF cos ip from the PI part [last line in (C2)]. 

Once the constants are chosen, the solution is uniquely fixed. The analysis of 
the solution u with the choice (C3) is well known and need not be discussed 
here. We could as well choose C\ = 1/R, Ci — in the zeroth order solution, 
and correspondingly the Pis would change, from which we can remove the CF 
in a similar manner as indicated above. 

As promised, let us consider the MKdS solution, first setting the trivial CFs 
to zero but including them in the Pis as they naturally appear. Then we get 

u = u + ui + u 2 (C4) 

where uq, u\ and ui are given by Eqs.(5), (9) and (10) respectively. For 7 7^ 0, 
this form of the solution is obviously different from that considered by Bodenner 
and Will because the last two Pis [u\ and 112) do contain terms proportional 
to trivial cost/?. [When 7 = 0, u consists of last three terms in (C2)]. The 
Bodenner- Will method has been applied to the asymptotically flat spacetime 
allowing the standard radial coordinate r — > 00. So they take u — > at tp — 
ir/2 + 5 and for small 6, sin 5 ~ S, cosS ~ 1. But u — > is not allowed in the 
present non-flat metric; instead one should tak^l u — u m in- Let us see what the 
Bodenner- Will method yields for the deflection angle. 

Ignoring the small terms proportional to S 2 , (-R7) 2 and (-R7) 3 in the equation 
resulting from (C4), we finally obtain 

«min = [4Mi?(8 - 37Ti? 7 ) - 8i? 2 (i? 7 + 26) 

-3M 2 {64i? 7 - US + 27r(4i?7(5 - 5)}] . (C5) 

This yields the value for 5 as 

_ 2(15M 2 tt + l&MR - 96M 2 i? 7 - 6MttR 2 1 - 4i? 3 7 ) - 16i? 3 w min 

_ 16i? 2 - 69M 2 + 24M 2 7ri?7 ' ' ' 

We straightforwardly expand 8 to find that 

R 3nM 423M 2 



2M 157rAf 2 
~R + 8R 2 7 
Y37T7 39 



V 2 16R 



2 4 32i? 

M 2 - r\ . (C7) 



2 The modification to r max of Case 2, Sec. IV by the presence of M is problematic, because 
the full equation B(r) = yields unwieldy roots. Hence we simply assume that there exists 
some r = r max or« = n min ~ 0. 
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Next, if we had excluded the CFs from PI (d la Bodenner and Will [15]), by 
suitably choosing constants we would obtain, under the same approximation, 



= — -t UMR(S - 371-.R7) - 8i? 2 (i?7 + 25) 
IqR 6 

-3M 2 {64i? 7 - 23(5 + 27r(4i? 7 (5 - 5)}] , (C8) 



which would yield 



r 2M 15ttM 2 



R 3ttM 453M 2 



2 4 32i? 



/3tt7 69 



^{{Y-m)"*-*} (C9) 

So in (C9), we find - 45 ^ instead of and -^f 3 ^ instead of 

— 39A f 6 ^ mi " , while other terms remain the same. That's all there is to it. How- 
ever, note that differences occur only in the third order of smallness, 7 ~^f = 

7.R x ^j-, which should not concern us as our purpose was to calculate effects 
only up to second order using the orbit equation to that order. Genuine third or- 
der effects would require the orbit equation in the third order in M. So we could 
as well delete those third order terms, none of which is genuine, from Eq.(21), 
(C7) and (C9) but we have nevertheless displayed them only to demarcate where 
the series ends. The remaining terms factored with u m j n are computed in the 
footnoteH The results of this Appendix therefore clearly demonstrate that dif- 
ferences in solution by inclusion or exclusion of trivial terms in PI (or choices of 
zcroth order sin <p, cos ip) lead to no differences in the final result for observable 
bending within the considered order. 

As stated in Sec. Ill, we wanted to faithfully adhere to the Rindler-Ishak 
[1] form of the solution (which they gave for 7 = 0). Also, we recover from 
Eq.(ll) the same form (sans lOsiny), now up to second order in M, as derived 
in Ishak et al [2]. They chose the definition of azimuthal angle (differing by ^ 
from the conventional one) that we also have maintained for the purpose of easy 
comparison with their expressions for - . Certainly, the final result for bending 
does not depend on this definition as shown in this Appendix, where the choice 
of the azimuthal angle has been conventional, as in Bodenner and Will [15]. The 
remarkable similarity of the terms up to second order in (C7) and (C9) above 



3 The exact value of r max (or « m i n ) from full B(r) = is rather messy. Hence, to simplify 
matters, let us assume that u m { n = — and that < R « ¥ =>■ kR << 7. With these 
assumptions, the first two terms factored with u m i n are at least of the third order of smallness, 
hence ignorable. The last one is — << 1 but its order of smallness depends on how small 

R is compared to An order of magnitude calculation assuming ? ~ 10 28 cm is in place. 



In the light-grazing-sun bending, numerology suggests that fc "^ ~ (iJq'J ' mcanm S third 
order of smallness, hence again ignorable. But for galactic clusters, say, Abell 2744, the values 
suggest that, while the first two terms are truly of the third order of smallness, the term 

— kRs ~ — 10 — 5 is in the first order. It is remarkable that the same first order effect follows 

1 

also from the Rindler-Ishak bending. 



1 3 
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with those in Eq.(21) shows that the bending derived via Rindlcr-Ishak method 
is indeed correct. 
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